Unusual spin-wave population in nickel after femtosecond laser pulse excitation 
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The spin-wave relaxation mechanisms after intense laser excitation in ferromagnetic nickel films are inves- 
tigated with all-optical pump-probe experiments. Uniform precession (Kittel mode), Damon-Eshbach surface 
modes and perpendicular standing spin waves can be identified by their dispersion co(H ext ). However, different 
to other ferromagnets co(H ext ) deviates from the expected behavior. Namely, a mode discontinuity is observed, 
that can be attributed to a non-linear process. Above a critical field the power spectrum reveals a redistribution 
of the energy within the spin-wave spectrum populated. 



I. INTRODUCTION 

Spin-wave generation and manipulation have already been 
demonstrated in many waysPJ They are of great importance 
in magnetism-based spin-wave logic and XNOR as well as 
NAND gates have recently been implemented. 2 To reduce the 
size of potential devices, effective spin-wave pumping mech- 
anisms are needed in ferromagnetic metallic films owning mi- 
cron to nanometer characteristic magnetic length scales. Here 
we present a pumping mechanism using fs-laser excitation 
that allows for a large precession angle (> 3°) after local ex- 
citation. On the other hand, the relevant length scales can be 
artificially selected by the design of magnonic crystals, i.e. 
periodically micro- and nano-structured magnetic materials, 
with only very selected dynamic eigen states. 

For large precession angles the equations of motion are in- 
trinsically non-linear. This can be used to pump energy into 
a certain spin-wave mode by non-linear interaction and be- 
comes evident for example in spin-wave nano-oscillatorsPS 
Given a spatial separation of two spin-torque oscillators and 
hence, an overlap in their spin-wave power when excited by 
a direct current, one finds a locking of frequencies if the fre- 
quency mismatch is below a critical value. Driven to an ex- 
treme, a steady transfer of energy from higher to lower en- 
ergy modes results in Bose-Einstein condensation observed at 
room temperature P 

In the following we focus on dipolar and exchange spin- 
wave modes in continuous films excited by intense laser 
pulses, where the nature of the excited spin-wave modes is 
tuned by the external magnetic field. The excitation by fem- 
tosecond laser pulses has major differences as compared to 
methods that work in thermal equilibrium. These are for ex- 
ample Brillouin light scattering (BLS)P conventional ferro- 
magnetic resonance (FMR), strip-line based techniques like 
vector network analyzer (VNAp or spatio-temporal Kerr ef- 
fect based FMR used to study spin- wave modes. 9 Our tech- 
nique works in space and time-resolved manner and the detec- 
tion process owns no ^-selectivity a priori. Yet, the strongest 
difference to reso nant te chniques is the broadband excitation 
in the spin system] 10 ! 11 ! by which all magnetic modes are ex- 
cited. However, only resonant modes contribute to the signal 
coherently and can be detected while others appear as non- 
coherent background. A further advantage is the contact-free, 
very local excitation and detection which simplifies the testing 
of various samples and structured media. 
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FIG. 1: Time-resolved measurement of the pump-induced change of 
the Kerr rotation (light gray) and the fitted background (dashed line) 
that is subtracted to yield the residual oscillation (black). The data 
set shown has been recorded at a nickel thickness of 80 nm under 
an external field of 140 mT. A large angle of precession as high as 
3° is observed. Inset: Schematic representation of spin-wave modes 
in a continuous ferromagnetic film of thickness d. Local amplitude 
and phase of the dynamic magnetization m(t) are illustrated. Shape 
anisotropy accounts for mainly in-plane alignment of the magnetiza- 
tion M = Ms +m(t). 



We make the very unusual observation, that below a criti- 
cal magnetic field two dominant spin-wave modes are excited, 
whereas above the critical field, they merge into one single 
mode. We shall also see, that intrinsic to the experimental 
method used, surface modes are preferentially excited if the 
sample thickness considerably exceeds the penetration depth 
of the laser field. 

In general, magnetic precessional motion can have dipo- 
lar and exchange spin-wave character and can be described 
as follows. Starting from the free magnetic energy density, 
the Landau-Lifshitz equation of motion for the magnetization 
vector M yields a theoretical dispersion describing the uni- 
form precession (also referred to as the Kittel mode) in the 
T3 
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FIG. 2: (color online). Spin-wave power spectra of thin nickel films with d = 80-180nm. Magnetic oscillations are observed after optical 
excitation (top plots) and normalization yields frequency branches (bottom plots), that can be attributed to uniform precession (white points), 
magnetostatic surface waves (blue points) as well as perpendicular standing spin waves (gray points). Solid lines represent theory according to 
Eqs. {T} to {3}. The vertical gray line at 80nm of nickel corresponds to the data set from |Fig.~T| vertical black lines for lOOnm denote Fourier 
spectra analyzed in |Fig. 4| Please note the mode discontinuity around /io^ext = 90 mT for d > 80 nm. 



where floM$ is the saturation magnetization (being = 659 mT 
for nickel), the only free parameter K z accounts for the ef- 
fective anisotropy in the out-of -plane direction, and H x is the 
component of the external field projected to the sample plane. 
Damon and Eshbach solved the equation of motion account- 
ing for Maxwell's equations in the magnetostatic limit and ad- 
equate boundary conditions yielding the dispersion for long 
wavelength surface modes, so called Damon-Eshbach (DE) 
modes: 



= H x (H x +M s ) + -f-[l- exp{-2k DE d)} . (2) 



Here, d is the film thickness and &de -L Ms denotes the wave 
vector of the magnetostatic surface wave. An analogous ap- 
proach in Ref. [[13 j yields dispersion characteristics taking 
into account exchange interaction as well as anisotropy. In 
thin magnetic layers one finds a quantization of the wave vec- 
tor perpendicular to the sample plane and in analogy to Eq. ([l} 
one can give a dispersion describing the exchange-dominated 
perpendicular standing spin waves (PSSW): 
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Therein, A is the exchange constant and £ pssw = n%d~ x is 
the quantized wave vector in out-of-plane (i.e. z)-direction at- 
tributed to a given order n of the PSSW. A schematic of the 
different magnetic modes is given in the inset of |Fig. 1| where 
the Kittel mode is not explicitly drawn, as both PSSW and DE 
modes result in the uniform precession in the limit fc pssw — > 
(i.e. n = 0) and &de — > 0, respectively. 



II. EXPERIMENT 

The samples consist of polycrystalline nickel films with 
thicknesses of up to 220 nm. The thickness dependence of 
the occurring modes is measured on a wedge-shaped sample 
to perform all measurements on one specimen. It was pre- 
pared by electron beam evaporation in ultra high vacuum un- 
der a base pressure of 5 x 10~ 10 mbar. A linearly moving 
shutter was used to produce a nickel wedge with thickness 
20 nm < d < 220 nm on a Si(100) substrate and in order to 
prevent oxidation, the ferromagnetic layer was capped with 
2 nm of copper. 

For the experiments we use an approach in the time domain, 
i.e. an all-optical pump-probe setup as previously described in 
Reference |fl4|. An ultrashort laser pulse with 60 fs duration 
and X c = 810 nm central wavelength excites the saturated sam- 
ple and at a variable time delay T a second laser pulse with the 
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FIG. 3: Spin-wave modes at fl()H ext = 50mT for varying thickness d. 
Below 40 nm only the uniform precession (white points, black line) 
is observed, which does not change in frequency for increasing d. 
Perpendicular standing spin waves with n = 1 can also be identified 
(gray) and fit well to the expected d~ 2 behaviorPS PSSW of second 
order are not observed (dashed gray line), instead Damon-Eshbach 
surface waves (blue) are excited. The blue dashed line represents the 
DE dispersion ^ with fcrjE = 1 um 1 an d reveals that the DE wave 
vector is not constant in our experiment (the solid blue line is a guide 
to the eye). 



same characteristics but 5% of the intensity probes the time- 
resolved magneto-optical Kerr effect (TRMOKE). 15 17 There- 
with, spin dynamics in the femtosecond regime can be de- 
tected, where electronic excitations lead to a demagnetization 
and subsequent relaxation into the initial state via scattering 
of high energy to low energy spin waves. 10 The processes that 
initially quench the magnetization in the first few picoseconds 
are subject to ongoing discussions, 18 but can be modeled quite 
successfully in some detail. 11 Here in focus are oscillations on 
a timescale from 30 ps to Ins. We shall see, that all modes 
corresponding to equations (jTJ to ([3]» can be optically excited 
in a thin ferromagnetic film. The responsible mechanism can 
be thought as an effective field pulse which stems from the 
heat-induced change of the sample's anisotropy and local ex- 
change field upon absorption of the pump pulse !^ Therefore, 
to promote magnetic precession the external field was tilted by 
30° out-of-plane thus giving rise to an angle between sample 
plane and the effective internal field Hj || Ms. For an external 
field of jUo^ext = 150mT, which is the maximum used in the 
experiments, we calculated a rotation of 7° of H, out-of-plane 
due to the Zeeman term in balance with the shape anisotropy 
The angle is rather small and thus, its influ- 



(see also Fig. 1 



ence can be neglected for simplicity. Hence, in the analytical 
expression of the spin-waves modes we assume //, and Ms to 
be in-plane. 



III. RESULTS AND DISCUSSION 

In order to obtain the precession frequencies from the time- 
dependent MOKE spectra, the incoherent background owing 
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FIG. 4: (color online). Gray and white points denote the Fourier 
power A^g as calculated from the time-resolved MOKE measure- 
ments at d = lOOnm, HoH exl = HOmT and 90mT, respectively. 
Solid and dashed lines represent Lorentz-Peaks which where fitted to 
the data. The expected linewidth in nickel is larger than 1 GHz and as 
a result, possible mode locking leads to a single peak at 1 lOmT if the 
overlap indicated by the hatched area is increased. Additionally, the 
mutual intensity is strongly increased. In the inset, the dependence 
of the Damon-Eshbach wave vector on the thickness is given, where 
the solid gray line is a guide to the eye (see text also). 



to phonons and magnons is subtracted. This procedure is 
demonstrated in |Fig. 1| with a single measurement recorded 
at /Xo/4xt = 140 mT and d = 80 nm. After subtraction of the 

Fourier transforma- 



background (black data points in Fig. 1 



tion of the M(t) curves yields peaks in the oscillations' power 
spectra. For a given thickness d, the external field is varied be- 
tween and 150 mT, resulting in a change of amplitude and 
frequency of the oscillation. The position of the correspond- 
ing peaks in the frequency domain is determined for each such 
set of measurements and plotted versus fioH ext (for details of 
the analysis refer to the supplementary data). 

Shown in |Fig. 2| are power spectra recorded on different 
nickel thicknesses, revealing up to three precessional modes 
of different origin. In the upper row the Fourier power spectra 
for a given thickness are plotted in a color map as calculated 
from the M(t) curves (the Fourier transformation of the data 
from Fig. 1 1 is marked with a vertical gray line). In the bottom 
row all spectra have been normalized by their respective max- 
imal FFT power to give a better overview of mode evolution. 
The data points represent the peak positions and included as 
solid lines are the fitted theoretical frequency dispersions ([T]i- 
<j3j of the various modes, which have been used to attribute 
the branches to the Kittel mode, the DE surface mode as well 
as the PSSW mode with n = 1. 

Analysis of the thickness dependent occurance of the dif- 



ferent magnetic modes yields the results depicted in Fig. 3 



For a given external field of 50 mT the observed frequencies 
have been extracted and the aforementioned identification of 
Kittel-, PSSW-, and DE-mode is confirmed. Both the Kittel 
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FIG. 5: (color online). Fourier power spectra for different pumping fluences F = 5-70mJcm~ 2 with d = 150nm. The identification of modes 
corresponds to |Fig. 2| and the color scale in the top row (not normalized plots) is the same for every spectrum to reveal the fluence-dependence 
of the precession amplitude. One observes a strongly increasing amplitude of the uniform precession with increasing fluence, whereas the 
bottom row (normalized data) shows a small shift of the critical field towards lower values, which will be detailed in |Fig. 6| 



mode (black) and the PSSW (gray) display the expected evo- 
lution with the tickness: While the frequency of uniform pre- 
cession Ok remains unchanged, the frequency of the PSSW 
reduces. Being proportional to d~ 2 it approaches »k in ac- 
cordance to literature ! 2 ™ 21 ! 



The DE modes (blue lines and points in |Fig. 3| > have been 
predominantly studied with Brillouin light scattering (BLS) 22 . 
We find from our experiment that they are first excited at a 
nickel thickness of 80 nm. There, they are only present at 
low external fields, but for thicker nickel layers these modes 
dominate the spin-wave dynamics over a wider field range up 
to fioHext = 90 mT. This is due to the excitation mechanism 
in combination with the relevant length scales in our experi- 
ment: sample thickness d and optical penetration depth A opt 
(i.e. penetration depth of the electrical field at optical frequen- 
cies which evaluates to A/(27l3(iV)) ~ 30nm, where N is 
the complex refractivity index). Calculations show that our 
experiment is sensitive to magnetic modes in the top 15nm 
of the ferromagnet, corresponding details will be published 
elsewhere.^ The magnetic excitation depth is governed by 
the penetration of the laser field and energy is deposited in 
form of an initial disarrangement of spins. 10 This energy will 
very effectively excite eigen modes with similar spatial pro- 
file and thus, by choosing a thickness d ^> A opt , the resulting 
strong asymmetry will make DE surface modes the favored 
relaxation channel. Consequently, the appearance of the DE 
modes is based on the fact that they own an amplitude pro- 
file with the maximum at the surface, decaying exponentially 
into the film. Herein lies a strong analogy to surface acoustic 
waves excited with laser pulses, which can be used to deter- 
mine elastic properties of surface layers!^ 



Assuming a propagation direction perpendicular to the 
magnetization (A:de -L Ms,H ext ) the dispersion |2]) can be used 
to determine the wave vector Icqe. In the inset of |Fig. 4 the 
respective quantitative analysis of the Damon-Eshbach modes 
is given. From the fits one obtains wave vectors in the in- 
verse micrometer range, corresponding to wave lengths of 
about 3um. We find that the product k^d entering Eq. Q 
is not constant but that the relation = (d — do)^ 1 holds 
with do = 66(4) nm (solid gray line in the inset of Fig. 4| . The 
divergence of &de around that value supports our earlier inter- 
pretation. An asymmetry between film thickness and optical 
excitation depth is needed to give rise to dipolar spin waves in 
contrast to exchange-dominated spin waves with rather high k. 
In other words, the parameter do can be considered an onset 
thickness above which the dipolar interaction comes into play. 
In general, our understanding of the selection of the DE wave 
vector is not fully developed. However, we can exclude, that 
the laser spot sizes of pump (60 pm) or probe (20 pm) enforce 
a selection, since the observed DE wave lengths are of the 
order of only very few microns and we do not observe any 
dependencies on the pump or probe beam diameter. 

As described, the DE mode dominates the spectra up to crit- 
ical fields as high as /Xo// C rit = 90 mT. Around that value, an 
explicit deviation from the theoretical dispersion arises (blue 
line and points in the bottom row of |Fig. 2) . A distinct modi- 
fication of the DE mode intensity takes place which - for suf- 
ficiently high fields - ultimately results in a peculiar charac- 
teristic evident in Fig. 2 Namely the DE mode merges into 



the Kittel mode. This leads to a very strong increase of the 
total Fourier power, i.e. precession amplitude, as seen in the 



non-normalized data (upper row of Fig. 2 1. The propagating 
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surface mode seems to lock to the frequency of the Kittel and 
PSSW modes for magnetic fields higher than IM)H cr i t . Above 
that value, the peaks can be attributed to the dispersion of the 
uniform precession by Eq. [T] The fits (white lines) yield val- 
ues for the effective anisotropy K z around 50kJm~ 3 showing 
no monotonous trend with the thickness. 

To discuss this unusual observation we want to remind the 
reader, that the highly asymmetric initial excitation profile un- 
favors the uniform precession as a magnetic eigen mode. One 
possibility to understand its observation is a locking of the DE 
mode to the k = homogenious precession originating from 
an energy transfer between the two dynamic modes. Only if 
the propagating surface waves would couple to the localized 
Kittel mode the amplitude of the magnetic precession could 
be that drastically increased. 

A reasonable starting point to discuss this idea is to draw 
parallelities to a model by Slavin and coworkers. They proved, 
that phase locking of spin-torque nano-oscillators can take 
placed In their publication, the authors use a non-linear set 
of equations of motion for two coupled oscillators and can ex- 
plain the experimental results presented in Refs. 14151 . Once 
the free running frequency mismatch between the two oscilla- 
tors is smaller than a threshold value A max (Eq. (9) in Ref. O), 
they phase lock to a mutual frequency. The mechanism is 
twofold: coupling can be either mediated by spin waves radi- 
ated into a common magnetic layer or a dipole field created 
by the oscillators. Either way, one of the crucial parameters is 
the spatial separation a of the two contacts, tuning the overlap 
of spin-wave power in real space (Fig. 1 in Ref. 0). 

Transferring these findings to our observation, we could 
speculate that an increased overlap of the magnetic modes in 
Fourier space leads to a sudden interaction such that locking 
takes place. Figure [4] illustrates this hypothesis: two Fourier 
spectra from measurements recorded on nickel with a thick- 
ness of lOOnm are shown. At jJ,oH ext = 90 mT clearly two 
precessional modes can be distinguished, whereas at /Xo#ext = 
HOmT only one smooth peak with nearly identical width is 
found. The two modes observed, namely uniform precession 
and DE surface waves, are separated in Fourier space by sev- 
eral GHz at low external fields. However, if jAoH ext is in- 
creased in the experiment, the separation reduces and at one 
point approaches the value of the linewidth, which implies 



an overlap in the frequency domain (hatched area in Fig. 4 1 



Above a threshold overlap the frequencies seem to lock and 
the two modes cannot be distinguished anymore. Moreover, a 
single smooth peak in the Fourier spectrum is observed. Con- 
cerning the evolution with increasing magnetic field, a step- 
like change of the DE frequency is apparent in the Fourier 
spectra in |Fig. 2| and hence, the term mode locking seems jus- 
tified. 

Equivalent to the damping time constant T a in the time do- 
main is the linewidth in the Fourier domain. From T a ap- 
parent in the time-resolved MOKE spectra one expects the 
width of the peak after Fourier transformation to be (nT a )~ l 
at 50% of maximum power (FWHM, full width half maxi- 
mum). For nickel the expected linewidth therefore is 1 . 1 GHz, 
using i a = 300ps, which has been extracted from the time- 
resolved spectra. Due to the strong damping in nickel (imply- 
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FIG. 6: (color online), (a) Threshold behavior of the Fourier 
power A^g jama °f tne Kittel mode. The external magnetic field 
can be used tune the observed characteristics between dipolar sur- 
face waves and uniform precession, (b) Determination of the critical 
magnetic field /io-f/ cl -;t for different fluences F — 5 and 60mJcm~ 2 , 
respectively. A^ e Kittel is plotted after normalization to the maximum 

of the respective Fourier spectrum for a given external field, A^ e max . 
The result can be regarded as a measure of mode population and 
a linear increase is found prior to saturation. The critical field can 
be extracted as depicted by the vertical green lines. In the inset, 
corresponding results for all fluences are included, giving a clear 
trend towards smaller /iof cl -; t for increasing F. Additionally, the Kit- 
tel amplitude as induced by different pump fluences is shown for 

W)# ex t = 150mT (g ra y point- 



ing a rather large line width) considerable overlap can occur. 
The interplay with H ext tuning the DE and Kittel frequency 
mismatch leads to possible phase locking (see |Fig. 4\ . 

A test for this hypothesis based on non-linearity is the influ- 
ence of the (optical) pumping power that excites the magneti- 
zation dynamics. Corresponding experiments were performed 
on a separate continuous nickel film which was 150nm thick 
and was prepared as described above. For the measurements, 
the pumping fluence was varied over one decade ranging from 
5 to 70mJcm~ 2 . The resulting Fourier power spectra are 
presented in |Fig. 5| and reveal that the nature of the mag- 
netic modes excited by the pump pulses does not significantly 
change. To give a better insight into the fluence dependence, 
all power spectra in the top row are identically color-scaled. 
This shows an increase of the spin-wave amplitude with the 
fluence (from left to right) as expected when considering the 
heat-induced disorder triggering the oscillation. 19 This in- 



crease is quantitatively given in the inset of Fig. 6 b), where 
the Kittel amplitude at /io^ext = 150mT has been plotted ver- 
sus the pumping fluence (gray points). In this respect, the 
absolute value of the amplitude is of great interest. It can be 
derived from the dynamic magnetization by calibration of the 
experiment for a given saturation magnetization: A compari- 
son of the magneto-optical signal at zero time delay with and 
without pump-induced demagnetization at F = 70mJcm~ 2 



6 



and /io^ext = 150mT results in a precessional amplitude of 
4.6% x floMs = 30mT, corresponding to an angle of 3°. 

A feature of |Fig. 5 1 that requires a closer look is the some- 
what changing critical field for the transition from surface to 
uniform characteristics. Figure[6|b) details the corresponding 
analysis and gives the respective results in the inset. Starting 
from the fit of the Kittel dispersion ([TJ to the data (white lines 
in the bottom row of Fig. 5 1 the amplitude of the Kittel mode 
for each magnetic field is extracted from the data. In the nor- 
malized case, this amplitude equals 1 if H ext > H cr \t and below 
the critical field a monotonous, in good approximation linear 
behavior is found. By a linear fit, the point at which the Kit- 
tel amplitude reaches 1, i.e. above which the Kittel mode is 
the dominant mode of precession can be determined. We de- 
fine this field value as the critical field // cr ; t (also marked by 
vertical lines in Fig. 6(b )). The above described procedure is 
explicitly given in Fig. 6 b) for the fluences 5 and 60mJcm -2 , 
respectively. In the inset the overall decrease of pk)H a \t for in- 
creasing F is shown. Together with the threshold-like onset 



of the uniform precession in Fig. 6 a) and the accompanying 



rapid increase of precession amplitude our initial hypothesis 
of a non-linear transfer of energy seems to be verified. 



IV. CONCLUSION 

In conclusion, we have identified the excitation of the 
uniform precession (Kittel mode) as well as magnetostatic 
surface waves (DE mode) in asymmetrically pumped nickel 



films. Below a critical magnetic field // cr j t the dipolar DE 
modes dominate the precession, whereas a step-like feature 
in the dispersion (£>Q^(H en ) indicates the onset of the pop- 
ulation of the k = mode above i/ cr j t (Fig. 2 1. The obser- 
vation of the DE surface waves for H ext < H CI i t can be well 
understood by considering the optical excitation mechanism. 
A concept was developed taking into account the amplitude 
profile of the DE modes as well as the asymmetric excitation 
profile stemming from the laser pump pulses. However, the 
understanding of the selection of relaxation channels in terms 
of the selection of the DE wave vector is a challenge still to be 
met. The transistion to uniform precession for H ext > H CI it was 
discussed in terms of a non-linear mechanism transferring en- 
ergy between the two eigen modes of the ferromagnet. Only 
a non-linear process allows to explain the dependence of the 
transition field // cr ; t on the pumping power. In this context the 
high energy of the pump pulses driving the system far out of 
equilibrium leading to large precession angles should be men- 
tioned. In total, the population of either the Kittel or the DE 
mode can be tuned by both the external field and the pump- 
ing fluence. Understanding the underlying mechanism might 
prove crucial on the way to spin-logic devices, especially with 
the interconversion of electrical, optical or magnetic signals in 
mind. 
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